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Abstract. Observations of active stars reveal highly non-uniform surface distributions of magnetic flux. Theoretical models 
considering magnetised stellar winds however often presume uniform surface magnetic fields, characterised by a single mag- 
netic field strength. The present work investigates the impact of non-uniform surface magnetic field distributions on the stellar 
mass and angular momentum loss rates as well as on th e effective Alfvenic radius of the wind. Assuming an axial symmetric 
and polytropic magnetised wind, the approach o f Web er & Davil J 19671) is extended to non-equatorial latitudes to quantify the 
impact of latitude-dependent magnetic field distributions over a large range of stellar rotation rates and thermal wind properties. 
Motivated by recent observational results, the analytically prescribed field patterns are dominated by magnetic flux concentra- 
tions at intermediate and high latitudes. The global stellar mass loss rates are found to be rather insensitive to non-uniformities 
of the surface magnetic field. Depending on the non-uniformity of the field distribution, the angular momentum loss rates devi- 
ate in contrast at all rotation rates between -60% and 10% from the Weber & Davisl -values. and the effective Alfvenic radii up 
to about ±25%. These large variations albeit equal amounts of total magnetic flux indicate that a classification of stellar surface 
magnetic fields through a single field strength is insufficient, and that their non-uniformity has to be taken into account. The 
consequences for applications involving magnetised stellar winds are discussed in view of the rotational evolution of solar-like 
stars and of the observational determination of their mass loss rates using the terminal velocity and ram pressure of the wind. 
For rapidly rotating stars the latitudinal variation of the wind ram pressure is found to exceed, depending on the actual field 
distribution on the stellar surface, over two orders of magnitude. The assumption of a spherical symmetric wind geometry may 
therefore lead to a significant over- or underestimation of the stellar mass loss rate. 
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1. Introduction 

Stars with hot coronae are subject to mass and angular mo- 
mentum losses due to stellar winds ( Parker 1958). In the pres- 
ence of magnetic fields the angular momentum (AM) loss is 
significantly enhanced, because the escaping ionised plasma 
is forced into co-rotation wit h the star out to d istances much 
larger than the stellar radius (Schatzman 1962). Further to its 
contribution to the specific AM of the plasma, the stress of 
bent magnetic field lines adds to the acceleration of the wind 
by magneto-centrifugal driving. In rapidly rotating stars this 
inherently latitude-dependent mechani sm is expected to domi- 
nate over the isotropic thermal driving (Michel 1969). 

Whereas early studies of stellar winds assumed uniform or 
dipolar surface magnetic fields, high-resolution observations of 
the solar atmosphere now indicate much more complex bound- 
ary conditions at the base of the wind. Doppler Imaging (DI) 
observations of active stars show non-uniform surface bright- 
ness distributions, which are ascribed to the presence of mag - 
netic flux in the form of dark spots ( Colli er Cameronl l2001 ). 
They show that on stars rotating more rapidly than the Sun 
magnetic flux is not only located in equatorial regions, but also 
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at intermediate and high latitudes (Strassmeier 2002, and ref- 
erences therein). Although starspots can a priori not uniquely 
be associated with wind loss regions, they nevertheless identify 
locations of possibly open magnetic field structures. Whereas 
spectroscopic DI enables the reconstruction of starspot maps, 
spectro-polarimetric Zeeman-DI directly confi rms the mag- 
netic origin of the dark surface features (Donati et alll997l) . In 
combination with field extrapolation techniques, the observed 
surface magnetic field distributions serve as boundary condi- 
tions for the reconstruction of coronal magnetic field topolo- 
gies and thus for the de termination of both closed and open 
field regions (e.g. Jardine et al. 2002). Recent results show that 
on rapidly rotating stars open magnetic flux is organised in belt- 
like structures main ly between intermediate and high latitudes 
jMcIvor et all2004l) . 

IWeber & Davisl ill 9671 hereafter WD) theoretically inves- 
tigated the radial structure and properties of the magnetised 
solar wind assuming that the wind structure determined in 
the equatorial plane is representative for the entire sphere, 
implying a monopolar, spherical symmetric field geometry. 
Parameter studies have been carried out in the framework of 
this model to investigate the dependence of the AM loss rate 
on the characteristic magnetic field str ength, the stellar rota- 
tion rate, and thermal wind properties ( Belcher & MacGregor 
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119761: Fya3 ll976). Mestel (1968) and IMestel & Spruid ill 9871) 
considered dipolar surface fields with closed magnetic struc- 
tures forming 'dead zones' at the e quator and radial o pen fields 
lines further away from the star; lOkamotol i 19741) extended 
this approach to different analytical poloidal field structures. A 
self-consistent determination of the poloidal field structure re- 
quires the solution of non-linear partial differential equations, 
which represents a formidable aspect of the overall wind prob- 
lem. Theoretical and numerical investigations including this 



of variables (Vlahakis & Tsinganos 


19981 iLimaetalJ 12001: 


Sautv et al. 


2002) or multi-dimensional MHD simulations 


(Sakurai 1985; 


Keppens & Goedbloedlll999. 


2000), respec- 



tively. They reveal at large distances from the star a rotation- 
dependent latitudinal deflection of the magnetic field lines to- 
ward the symmetry axis as a result of the collimating e ffect 
of the azimuthal field component (e. g . iHe vvaerts & Nor man 
Il989t iTsineanos & Bogovalovl Eoool lOkamotol l2000l) . The 
self-consistent treatments are mainly focused on the (asymp- 
totic) wind structure without going into details about mass and 
AM losses for different stellar parameters. Owing to the non- 
linearity of the influence of the poloidal field component on the 
overall wind structure, the differences between self-consistent 
methods and simplified approximations are difficult to predict 
and have to be verified through direct comparisons. Therefore, 
the results of WD-like approaches should a priori be regarded 
as estimates within the framework of the applied simplifica- 
tions. 

The appro a ches of IWeber & Davisl i 19671) and 
IMestel & Spruid Jl987t> are widely used for the determi- 
nation of A M loss rates in the cour s e of the rotational evolution 
of s tars |Endal & Sofial 1 198 it ICoUier Cameron & Tiankel 
1994 1 O Il999l). I n cont rast, ISolanki et alJ i 19971) and 
Holzwarth & JardinJ J2005I) used specific modifications of 
the WD model to investigate the impact of non-uniform 
surface fields. To determine m ass loss rate s of sol ar-like stars, 
IWood & Linskvl J 19981) and IWood et alJ d2002h devised a 
method using the line-of-sight absorption of the Lya emission 
by surrounding astrospheres, which are formed by the collision 
between the (magnetised) wind from the star and the local 
interstellar medium. The size of the astrosphere as well as the 
amount of absorption depend on the ram pressure of the wind, 
which is a function of the stellar mass loss and wind velocity. 
For the determination of the total stellar mass loss rate they 
assume a constant (i.e. spherical symmetric and independent 
of the rotation rate and latitude) wind velocity. Based on their 
results, they formulate an empirical relationship between mass 
loss rates and coronal X-ray fluxes of solar-like stars. 

Motivated by the progress in the determination of stellar 
surface features and mass loss rates, the present work applies 
the WD wind model to non-equatorial latitudes to investigate 
the influence of latitude-dependent surface magnetic fields on 
the total mass and AM loss rates of cool stars. In contrast 
to self-consistent treatments of the magnetised wind problem, 
which are focused on detailed properties of the wind structure, 
the main objective here is the global qualification and quantifi- 
cation of the influence of observed surface magnetic field dis- 
tributions in relation to the originally spherical symmetric WD 



method. The applied wind model represents an efficient way to 
include the aspect of non-uniform surface magnetic fields in the 
rotational evolution of solar-like stars and in the observational 
determination of their mass loss rates. 



2. Model assumptions 

The stellar mass and AM loss rates through latitude-dependent 
magnetised winds are determined for a solar-like star by ap- 
plying the WD approach to open magnetic field lines at all lat- 
itudes separately (cf. Appendix 1X1. The poloidal components 
of the magnetic field and the flow velocity are assumed to 
be radial, neglecting the trans-field component of the (station- 
ary) equation of motion perpendicular to the magnetic field. 
Retaining the full azimuthal wind dependency, the open field 
lines form spirals on coni with constant opening angles, whose 
tips are located in the stellar centre. 

The polytropic wind is assumed to be symmetric with re- 
spect to both the rotation axis and the equator, enabling a limi- 
tation of the analyses to the range of co-latitudes < 9 < n/2. 
The wind solutions are determined through the stellar rota- 
tion rate and the magnetic and thermal wind conditions pre- 
scribed at a reference level, ro = 1.1 Ro, close to the stel- 
lar surface. The uniform rotation rates cover the range Q = 
10~ 6 - 5.4 ■ 10 _4 s _1 ~ 0.3 - 180 £2 Q , that is rotation periods 
73 d - 3 h; the upper limit of the rotation rates corresponds to 
the case when the co-rotation radius, Eq. JA.12> . in the equato- 
rial plane is equal to the reference radius. 



2.1 . Thermal and magnetic wind properties 

The wind temperature and (particle) density at the reference 
level are approximated through solar-like values, that is Tq = 
2 ■ 10 6 K and no = 10 8 cirT 3 , respectively. With a mean atomic 
weight = // a 0.6, these values imply a gas density of 
po — 10~ 16 g/cm 3 and a gas pressure of /?o - 3 ■ 10~ 2 dyn/cm 2 . 
The entropy change of the wind with increasing distance from 
the star is quantified through the polytropic index, T = 1.15. 
These values fall in the range of values of simil ar studies (e.g., 
r * 1.2 and T = 2.7 • 10 6 K. IWeber & Davislll967l ISakurail 
Il985t T = 1.13 and T = 1.5 ■ 10 6 K. lKeppensetalJll995l) . 
and have been chosen to ensure that even plasma emanating at 
high latitudes, where the magneto-centrifugal driving is inher- 
ently less efficient, escapes from the star owing to a sufficiently 
strong thermal driving. The thermal boundary conditions are 
taken to be independent of both latitude and rotation rate. 

The magnetic wind properties are determined through the 
radial magnetic field component at the reference level. For the 
uniform Constant Field (CO) distributions the field strengths 
are taken to be B r $ = 3,30, and 300 G, corresponding to an 
(unsigned) open magnetic flux between cD =s 2 ■ 10 23 ~ 25 Mx. 
The lower limit is of the order of the average solar mag- 
netic field strength, whereas the upper limit is in accor- 
dance with spectro-polarimetric observations of rapidly rotat- 
ing stars, which indicate the existence of field strengths up to 
two orders of magnitude larger than in the case of the Sun 
dDonati & Collier CamerorJ 19971) . 



V Holzwarth: Latitude-dependent magnetised winds 



3 



AB 



PS 




Fig. 1. Polar Spot (PS) and Latitudinal Belt (LB) distributions 
of the radial magnetic field component, B r fl{9), at the reference 
level. 



In the case of non-uniform magnetic fields, latitude- 
dependent distributions are superposed on a constant back- 
ground field of B < = 3 G (Fig. 0. The peak field strength, 
B > — B K + AB, is chosen so that the surface averaged mag- 
netic field strength, 



tt/2 



B 



r0 = J B rfi (0) sin dO, 
o 



(1) 



of each distribution is 30 G, to allow for a comparison with 
the intermediate case of the Constant Field distributions above. 
Using the analytical form 



B rfi (0) - B < + AB ■ cos" (0 - 6> ) 



(2) 



the latitude-dependent magnetic field distributions considered 
here are the 



8,0o = 0, and AB 



1. Polar Spot model (PS), with n 
243 G, and the 

2. Latitudinal Belt model (LB), with n = 16, O = 45° 
Afi = 63.8G. 



and 



The first distribution is motivated by observati onally deter- 
mined surface brightness maps of active stars JStrassmeieJ 
120021) . which frequently show spot concentrations at high 
latitudes. The second one is initiated through recent results 
of coronal magnetic field extrapolations based on Zeeman- 
Doppler Imaging observations of rapid rotators, which reveal 
a gathering of open magnetic flux also at intermediate latitudes 
jMcIvoretalJ2n04l). 

2.2. Stellar mass and angular momentum loss rates 

The contribution of a plasma stream to the stellar mass loss rate 
is given by 



dM = F wider = 2ni r A p A v r ^ A sin d6 . 



(3) 



with dcr = sin d8 dip being the solid angle occupied by the 
escaping plasma flow. The mass flux per solid angle, Fm = 
pv r r 2 = p A v r ,Ar A , ' s constant along individual field lines but 
latitude-dependent, since the plasma density, pa, and radial 



flow velocity, v r , A , at the Alfvenic point, r A , are functions of 
the co-latitude. The total stellar mass loss rate is given by 



M — Mwd fN ( P A(—V^d6. 

J W/ \PKl\%Al 



The WD mass loss rate, 
M W d = 47rr A p A v,-,A , 



(4) 



(5) 



is derived from the assumption that the density, pA, and radial 
flow velocity, v r _A, at the Alfvenic point, ?a, in the equato- 
rial plane are representative for the wind structure at all lat- 
itudes, implying r A (0) = r A ,p A (0) = Pa, and v r , A (0) = v r ,A. 
All WD values are calculated using the surface averaged field 
strength, Eq. ([Q; the determination of wind solutions is out- 
lined in AppendixlAl 

The AM momentum loss per latitude is 



dj — L Folder — 2n£lr A p\v r ^ sin 3 d6 , 



(6) 



where L — Q (r A sin 0) is the specific AM along an individual 
field line [cf. Eqs. $K3l & Q ], The total AM loss rate of the 
star is then given by 



n/2 



1 - 



with the respective WD value 



87T a _ 2 2 • 

/wd = — i2r A p A v r A = -Or A M W D 



(7) 



(8) 



Latitude-dependent boundary conditions at the base of the stel- 
lar wind result in quantities (r A ,p A , v nA ), which are functions 
of the co-latitude and a priori different from (r A ,p A , v r jd in the 
equatorial plane. Even boundary conditions constant at all lat- 
itudes cause a deviation of the wind structure from spherical 
symmetry, since the contribution of the magneto-centrifugal 
driving to the overall wind acceleration along non-equatorial 
field lines is smaller than in the equatorial plane. 

Studies considering the rotational evolution of stars of- 
ten determine AM loss rates not through explicit boundary 
conditions, but by using Eq. © in combinati on with either 
prescribed or emp irical mass loss rates (e.g., lKawalerl ll988: 
IWood et aD l2002). deriving the required Alfvenic radius by 
more qualitative arguments and/or relationships. The influence 
of non-uniform field distributions on this approach is quantified 
through the effective Alfvenic radius, 



>"A,eff = r A 



(Ml 



WD 



WD 



1/2 



(9) 



which describes the functional relation between the latitude- 
integrated mass and AM loss rates in comparison with the WD 
results. 
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Fig. 2. Reference values determined following the approach 
of IWeber & Davisl ll 19671) . with B rfi = 3 G (solid), 30 G (long 
dashed), and 300 G {short dashed). Panel a: Mass loss rate, 
Mwd; for the two larger field strengths the mass loss rates are 
quite similar and hardly distinguishable. Panel b: Angular mo- 
mentum loss rate, Jwd', the thin line indicates the slope of a 
function linear in the stellar rotation rate, SI. Panel c: Specific 
angular momentum, Lwd- Panel d: Alfvenic radius, twd; the 
thin line indicates the slope of a function oc ST 1 ! 1 . 

3. Results 

3.1. Reference cases 

The mass loss rates, Mwd, AM loss rates, 7wd, and Alfvenic 
radii, f\, determin ed in the equa torial plane following the ap- 
proach of IWeber & Davisl Jl967l) . are used as reference values. 
The stellar mass loss rate is controlled by the wind structure 
in the subsonic flow regime close to the star. Its value depends 
significantly on the stellar rotation rate but hardly on the mag- 
netic field strength (Fig.|2k). The increase of the AM loss rate 
with the rotation rate (Fig. |3>) is due to the strong increase 
of the mass loss rate, which outbalances minor variations of 
the specific AM (Fig. |2j;). Approximations based on rotation- 
independent magnetic field strengths imply a linear increase 
of the AM loss rate, that is j oc SI" with a =s 1, whereas the 
functional relations here are found to be sub-linear. For rotation 
rates SI > 2 ■ 10~ 5 s _1 ^ 7Qq a curve fit yields a 0.8. Below 
this rotation rate a =* 0.8 (3 G), 0.45 (30 G) and 0.4 (300 G), 
respectively. Due to its explicit linear dependence on the stellar 
rotation rate, the specific AM increases at slow rotation rates, 



Table 1. Refe rence valu es determined following the approach 
of IWeber & Davisl dl967h . for B n0 = 30 G. The solar rotation 
rate is taken to be Qq = 2.9 • 10~ 6 s _1 . 
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Fig. 3. Relative deviations, 6M, of the latitude-integrated mass 
loss rates from the WD reference values. Constant Field [solid); 
Polar Spot {long dashed); Latitudinal Belt (short dashed) dis- 
tribution with B r o = 30 G. The thin grey shaded region indi- 
cates the range of values for Constant Field distributions with 
field strengths between 3 - 300 G. Arrows mark rotation rates 
of the latitude-resolved cases shown in Fig.0] 

where the decline of the square of the Alfvenic radius is sub- 
linear (Fig. EJj). At higher rotation rates this decline becomes 
super-linear and the specific AM consequently smaller, entail- 
ing a maximum at intermediate rotation rates. 

To verify the impact of non-uniform field distributions with 
respect to the spherical symmetric WD approach, in the follow- 
ing the results are mainly given in terms of relative deviations, 
6F — (F — Fwd) /Fwd, where F is the mass loss rate, Eq. @, 
the AM loss rate, Eq. @, or the (effective) Alfvenic radius, 
Eq. 0, respectively. Further quantities have been used to nor- 
malise the latitudinal profiles in Figs.^lU and|SJ their values 
are given in Tabled 

3.2. Mass loss rate 

With increasing rotation rate the latitude-integrated mass loss 
rates fall short of the WD value, in the domain of rapid rotators 
up to about 35% (Fig. [3}. The major mass loss originates from 
the equatorial region, where the magneto-centrifugal driving 
of the wind is most efficient. However, the mass loss rate per 
latitude, dMjdQ, decreases significantly with increasing lati- 
tude (relative to the WD values), depending on the rotation rate 
about two to three orders of magnitude (Fig. |4ji. This strong 
decline is caused by the decreasing mass flux in combination 
with the latitudinal weighting function oc sin 8, which takes the 
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Fig. 4. Left quadrant: Relative mass loss rate per latitude, 
lg ((dM/d8)/(dM/d6)y/r,\ Right quadrant: Relative mass 
flux, lg (F M /F M WD ). Polar Spot distribution with £2 = 2.9 • 
1(T 6 s- 1 (solid); 2.2- 1(T 5 s" 1 (long dashed); 1.4- 1(T 4 s" 1 (short 
dashed); the respective latitudinal profiles of the Latitudinal 
Belt and Constant Field distributions are qualitatively very sim- 
ilar and therefore not shown. Note that with r^po = const, the 
relation Fm <* v r ,o holds. 



smaller surface fraction at higher latitudes into account. In the 
regime of slow stellar rotation the wind acceleration is caused 
by the latitude-independent thermal driving, so that the mass 
flux is virtually constant at all latitudes. In this case the latitudi- 
nal variation of dM/d9 is dominated by the weighting function. 
For larger rotation rates the latitudinal decrease of the mass flux 
strongly enhances the drop of the total mass loss rate. Since the 
WD approach implies that the equatorial wind structure is rep- 
resentative for all latitudes, the generalisation of the equatorial 
mass flux to higher latitudes results in an overestimation of the 
WD mass loss rate in comparison with the latitude-integrated 
values. 

The deviations of the mass loss rates subject to latitude- 
dependent field distributions are larger than in the case of con- 
stant surface fields, but retain the same functional behaviour, 
with largest deviations occurring at high rotation rates. The 
Polar Spot yields mass loss rates up to ~ 10% smaller than 
the Constant Field distributions, whereas the Latitudinal Belt 
falls short only by a few percent. High-latitude magnetic flux 
concentrations like a Polar Spot imply smaller field strengths at 
equatorial regions and consequently a drop of the mass flux and 
mass loss rate per latitude. The Latitudinal Belt possesses at in- 
termediate latitudes in contrast sufficiently large field strengths 
and mass fluxes to sustain mass loss rates similar to the case 
of a Constant Field. Since non-equatorial latitudes contribute 
less to the overall mass loss rate, this localised mass flux sur- 
plus can however not balance the reduction of mass loss at the 
equator. 

3.3. Angular momentum loss rate 

The AM loss rates depend significantly on the stellar rotation 
rate and magnetic field strength, but even stronger on the lat- 
itudinal field distribution (Fig.|5Jl. In the case of constant sur- 
face fields, the deviations of the latitude-integrated AM loss 
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Fig. 5. Relative deviations, 6 J, of the latitude-integrated an- 
gular momentum loss rates from the WD reference values. 
Constant Field (solid); Polar Spot (long dashed); Latitudinal 
Belt (short dashed) distribution with B r ,o = 30 G. The grey 
shaded region indicates the range of values for Constant Field 
distributions with field strengths between 3 - 300 G. Arrows 
mark rotation rates of the latitude-resolved cases shown in 
Fig. El 

rates from the WD values are at small rotation rates moderate, 
up to about 10%, and field strength-dependent. In contrast, at 
large rotation rates the values become only marginal and inde- 
pendent of the field strength. In the case of non-uniform field 
distributions, the integrated AM loss rates can be significantly 
different despite of equal total amounts of open magnetic flux. 
In the case of a rapidly rotating star with a Polar Spot the total 
AM loss rate is about 60% smaller than the WD value. In the 
case of a Latitudinal Belt distribution it is about 10% larger. 
In contrast to the mass loss rate, where the deviations increase 
from marginal values at small rotation rates up to a maximum 
in the regime of fast rotation, the deviations of the AM loss rate 
are significant over the entire range of rotation rates. The func- 
tional dependence on the stellar rotation is similar for all field 
distributions. 

The functional behaviour of the AM loss rate is caused by 
its combined dependence on the mass loss rate (per latitude) 
and the specific AM, the latter being very susceptible to field 
strength variations (Fig.|6j. The specific AM is a latitudinally 
weighted function oc sin 2 6, due to its dependence on the axial 
distance of the Alfvenic point. The Alfvenic radius itself in- 
creases with latitude and decreases with the rotation rate (cf. 
Fig. |2ji, which reflects the changing efficiency of the magneto- 
centrifugal driving. In the case of uniform magnetic fields, the 
foreshortening of the geometric lever arm toward high lati- 
tudes is therefore countered by the concomitant increase of 
the Alfvenic radius, and a dependence of the specific AM on 
the rotation rate opposite to the one of the mass loss rate. The 
combination of the latter two quantities results in a lateral com- 
pensation and therefore in similar t/y/iifl-profiles (Fig. |6j CO). 
Consequently, the WD approach underestimates the specific 
AM of non-equatorial regions by assigning the equatorial value 
of the Alfvenic radius to all latitudes. This underestimation is 
partly balanced by the overestimation of the mass loss rate, 
which eventually entails moderate deviations of the latitude- 
integrated AM loss rates from the WD values. 
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Fig. 6. Left quadrants: Relative angular momentum loss 
rate per latitude, lg [(dj/d9)/(dj/d6)wD\ Right quadrants: 
Relative specific angular momentum, lg[L/LwD]- Constant 
Field (CO); Polar Spot (PS); Latitudinal Belt (LB) distribution 
with B rfi = 30 G and Q = 2.9 ■ 1(T 6 s" 1 (solid); 2.2 ■ 1(T 5 s" 1 
(long dashed); 1 .4 • 1CL 4 S _1 (short dashed). The dashed-dotted 
lines show the weighting function oc sin 2 9; since L oc r 2 sin 2 6, 
the differences between this line and the solid/dashed lines are 
caused by the latitudinal variation of the Alfvenic radius. 

A Polar Spot distribution yields the largest differences in 
the latitudinal profile of the specific AM. Although the mag- 
netic field near the pole is very strong, its contribution to the 
AM loss rate (per latitude) is only average, since at high lati- 
tudes the mass loss rates are very small. The reduction of the 
specific AM at low latitudes causes in contrast a much larger 
deficit, so that the total AM loss rate is far below the WD value. 
The Latitudinal Belt distribution also shows a drop of the spe- 
cific AM around the equator, but since its low-latitude field 
strengths are larger than in the Polar Spot distribution the differ- 
ences are in total rather small. The concentration of magnetic 
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Fig. 7. Relative deviations, 5r ae (f, of the effective Alfvenic 
radii from the WD reference values. Constant Field (solid); 
Polar Spot (long dashed); Latitudinal Belt (long dashed) dis- 
tribution with B = 30 G. The grey shaded region indicates 
the range of values for Constant Field distributions with field 
strengths between 3 - 300 G. 

flux in a Latitudinal Belt implies high magnetic field strengths 
at intermediate latitudes. The large specific AM yields in com- 
bination with the substantial mass loss rates AM loss rates (per 
latitude), which are even larger than in the reference case. 

3.4. Effective Alfvenic radius 

For the cases considered here, the deviations of the effective 
Alfvenic radius from the WD values are within about +25% 
(Fig-0- For a Constant Field or Latitudinal Belt distribution the 
effective Alfvenic radius is larger than the WD value; at small 
rotation rates the difference is only marginal, but in the case 
of rapidly rotating stars it reaches a maximum of ~ 20%. The 
Polar Spot distribution, in contrast, yields effective Alfvenic 
radii which are for all rotation rates nearly 25% smaller than 
the reference values. 

3.5. Terminal wind velocity and ram pressure 

The ram pressure, p w , exerted by a magnetised stellar wind far 
away from the star is approximately proportional to the prod- 
uct of its mass flux (cf. Fig.|4} and terminal wind velocity, v M . 
The latitudinal profile of the latter reflects the strong depen- 
dence of the magneto-centrifugal driving on the magnetic field 
distribution and on the stellar rotation rate (Fig.|8|i. For small 
field strengths and/or rotation rates this driving mechanism is 
rather inefficient and the terminal wind velocity thus small. In 
the case of slowly rotating stars, the latitudinal variation of the 
ram pressure is found to be within a factor of ~ 3 compara- 
ble with the WD value, /? w ,wd k A^wdVm.wd- At high rotation 
rates, however, the latitudinal variation depends significantly 
on the latitude-dependent magneto-centrifugal driving of the 
wind, and consequently on the actual magnetic field distribu- 
tion. In case of a Polar Spot the wind acceleration at high lat- 
itudes is supported by strong magnetic fields. There the large 
terminal velocities make up for smaller mass fluxes, so that the 
resulting latitudinal variation of the ram pressure is about one 
order of magnitude. In case of the Latitudinal Belt distribution, 
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Fig. 8. Left quadrants: Relative terminal wind velocity, 
lg (voo/v^wd)- Right quadrants: Relative wind ram pres- 
sure, lg (Pw/Pw/wd)- Constant Field (CO); Polar Spot (PS); 
Latitudinal Belt (LB) distribution with B rfl = 30 G and Q = 
2.9- 10-V 1 (solid); 2.2 • 1(T 5 s" 1 (long dashed); 1.4- lO^s" 1 
(short dashed). 



for which the magnetic flux is concentrated at intermediate lat- 
itudes, the terminal wind velocities at high latitudes become 
very small and the difference of the ram pressure between in- 
termediate and polar latitudes very large, here over two orders 
of magnitude. 

4. Dependence on thermal wind properties 

The winds considered here are accelerated through thermal and 
magneto-centrifugal driving. The latter mechanism is in con- 
trast to the isotropic thermal driving inherently rotation- and 
latitude-dependent. The fractional contribution of each driv- 
ing mechanism to the overall wind acceleration determines 
the emphasise of its latitude-dependent character, which con- 



sequently changes with the stellar rotation rate. A priori, the 
thermal driving complies better with the basic WD assump- 
tion of a spherical symmetric wind geometry than the latitude- 
dependent magneto-centrifugal driving. To quantify the im- 
pact of the thermal wind properties on the previous results, 
each one of the thermal parameters is successively changed, 
while the others are kept constant; the parameter ranges are 
T {) = (1.7-2.5)- 10 6 K;« = 10~ 7 -10~ 9 ctrT 3 ; F = 1.05-1.175. 

The mass loss rates of rapidly rotating stars, whose out- 
flows are dominated by the magneto-centrifugal driving, are 
generally less susceptible to changes of the thermal wind prop- 
erties than those of slow rotators. Higher temperatures and/or 
stronger heating (i.e. smaller polytropic indices) of the wind 
shift the emphasise of its acceleration from the magneto- 
centrifugal to the thermal driving (Fig. [5] top row). The transi- 
tion between the regimes of thermally controlled and rotation- 
ally controlled winds thus occurs at somewhat higher rotation 
rates. Since the thermal driving is taken to be isotropic, the inte- 
grated mass loss rate becomes in turn less susceptible to latitu- 
dinal variations of the magnetic field, so that the relative devia- 
tions, 6M, are smaller. For cooler wind temperatures the effect 
is vice versa; the mass loss rates are up to 40% smaller than 
the WD reference values, with considerable deviations even at 
small rotation rates. Tenuous winds with lower densities result 
in smaller deviations of the mass loss rate than more substan- 
tial outflows. Thereby, magnetic field concentrations at high 
latitudes entail a larger susceptibility of the mass loss rate to 
variations of the wind density than field concentrations at lower 
latitudes. 

The impact of the thermal wind properties on the AM loss 
rate is rather complex and cannot easily be described in sim- 
ple terms (Fig. [9] middle row). The specific AM along an open 
field line depends on the axial distance of the location where 
the Alfvenic Mach number (oc v yfp/B) is unity. Therefore, the 
smaller the magnetic field strength, or the larger the wind den- 
sity and/or flow velocity, the closer to the star the Alfvenic 
point will be located. For higher wind temperatures the spe- 
cific AM is therefore smaller. But since this decrease is out- 
balanced by the concomitant increase of the mass loss rate, 
the resulting AM loss rates are usually larger than for cooler 
winds. For the thermal wind parameters considered here, the 
(^./-deviations have magnitudes roughly similar to the standard 
cases (thick solid lines in Fig.|5}, although the individual func- 
tional behaviour can be quite different. In particular at small 
rotation rates the deviations are susceptible to the thermal wind 
properties, showing a rather non-linear behaviour. 

Owing to its definition, Eq. 0, the effective Alfvenic ra- 
dius reflects the functional dependence of both the mass and 
AM loss rates. Since the deviation of the mass loss rate is typ- 
ically larger than the deviation of the AM loss rate, the func- 
tional dependence of (5rA, e ff is less complex than 6J; an excep- 
tion is the case of the Polar Spot distribution, for which the im- 
pact on the AM loss rate is very large. The dependence on the 
wind density is generally rather small, whereas cooler (hotter) 
winds are found to entail larger (smaller) deviations from the 
WD values. The same is respectively true for the dependence 
on the wind heating. 




Fig. 9. Relative deviations of the mass loss rate (top), the AM loss rate (middle), and the effective Alfvenic radius (bottom) for 
the Constant Field (CO); Polar Spot (PS); Latitudinal Belt (LB) distribution. The thermal properties at the base of the wind are 
r = 1.15 and To = 2-10 6 K = T lef ,n = 10 s cnr 3 = n ref (solid); T = 2.5-10 6 K,n = n ief (long dashed); T = 1.7-10 6 K,« = n lef 
(short dashed); To — T re f,rio = 10 9 cirT 3 (dashed-dotted); To = T re {,no = 10 7 cirT 3 (dashed-triple dotted). The grey shaded 
regions indicate the range of values for polytropic indices T — 1.05 - 1.175 and T le f, n, e f; the impact of a lower polytropic index 
is thereby comparable with a higher wind temperature. The (surface averaged) magnetic field strength is in all cases B r ,o = 30 G. 



5. Discussion 

The non-uniform field distributions considered here are moti- 
vated through observations of rapidly rotating stars, which fre- 
quently indicate the occurrence of magnetic fields up to polar 
latitudes. Theoretical models explain these high-latitude fea- 
tures in terms of the pre-emptive poleward deflection of mag- 
netic flux by the Coriolis force (Schiissl er & Solankil ll992). 
and/o r its post-eruptive p oleward transport by meridional mo- 
tions ( Sch riiver & T itle 2001 ). The approach of Weber & Davis 
(1967, WD) in contrast implies a split monopole (i.e. spher- 
ically symmetric) geometry of the magnetic field and wind 
structure, which results in an overestimation of the mass flux 
and an underestimation of the specific AM at high latitudes. 
The latitude-integrated mass and AM loss rates subject to non- 
uniform surface fields show consequently large deviations from 
the WD approach. 

The deviations of the mass loss rate are largest (up to 
30 - 35%) for rapidly rotating stars, since in this regime the 
dominance of the latitude-dependent magneto-centrifugal driv- 
ing causes large asymmetries in the mass flux between the 
equator and the pole. The dependence of this quantity on the 
field strength is small, so that a redistribution of magnetic flux 
from equatorial to intermediate or high latitudes hardly changes 
its latitudinal pattern. The mutual deviations of the mass loss 



rates between different non-uniform flux distributions are here 
< 10%. 

Since the functional dependencies of the mass flux and of 
the specific AM on the rotation rate follow opposite trends, 
their combined effects partly balance, which makes the devia- 
tion of the AM loss rate from the WD value nearly rotation- 
independent. Non-uniform fields, on the other side, impose 
considerable variations on the latitudinal distribution of the 
specific AM, which are inherited by the AM loss rate caus- 
ing deviations between -60% and 10%. Typically, the stronger 
the concentration of magnetic flux at higher latitudes, the 
smaller the AM loss rate. However, for flux concentrations in 
a Latitudinal Belt at intermediate latitudes the AM loss rate 
is found to be larger than the WD value, because the higher 
magnetic field strengths outbalance the foreshortening of the 
geometric lever arm. 

The effective Alfvenic radius reflects the individual proper- 
ties of the mass and AM loss rates; its deviation form the WD 
value are here found to be within about ±25%. Owing to the 
quadratic dependence on the Alfvenic radius, the AM loss rate 
can thus deviate up to +50% from the WD values, if it is de- 
termined in terms of an empirical mass loss rate. The largest 
uncertainties in the AM loss rate caused by disregarding the 
actual non-uniformity of surface magnetic fields are likely to 
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occur in the case of rapid stellar rotation, when magnetic flux 
is gathering at high latitudes. 

Both higher wind temperatures and densities are found to 
increase the mass loss ra te, in agreement with previous results 
jMestelll968l:lYehll976l) . With increasing wind temperature (or 
heating) the isotropic thermal driving becomes stronger, and 
the importance of the magneto-centrifugal driving and of non- 
uniform field distributions weaker. The transition from ther- 
mally to rotationally controlled winds is thus shifted to higher 
stellar rotation rates. Vice versa, the mass loss rates of cool 
winds are more susceptible to non-uniformities of the sur- 
face magnetic field. The AM loss rates as well as the effec- 
tive Alfvenic radii show a more complex dependence on the 
thermal wind properties. The order of magnitude of their de- 
viations from the respective WD values remains roughly unaf- 
fected, but the functional behaviours change in complex ways. 
In accordance with the solar paradigm, high coronal tempera- 
tures are ascribed to enhanced magnetic activ ity, which is fo und 
to in crease with the stellar rotation rate dMontesinos et alJ 
l200lh . Consequently, one has to assume that in rapidly rotating 
stars the structuring of stellar winds arising from the latitude- 
dependent magneto-centrifugal driving and from non-uniform 
surface magnetic fields is diluted by high wind temperatures, 
which in turn fortify the isotropic thermal driving. 

In this work, the range of magnetic field strengths of 
the non-uniform field distributions is bracketed through the 
low- and high-field strength cases of uniform surface fields. 
However, as mentioned above, their individual loss rates are 
significantly different. According to this, neither the surface 
averaged, nor the maximum, nor the minimum field strength 
enable an accurate characterisation of stellar magnetic fields. 
In contrast, the AM loss rates of different non-uniform field 
distributions can cover a large range of values, although their 
total magnetic flux through the stellar surface is the same. This 
has important consequences for studies considering the rota- 
tional history of stars and the distribution of their rotation rates. 
By investigating the impact of non-unifor m surface mag netic 
fields o n the rotational evolution of stars, Holzwarth & Jardine 
(2005) show that concentrations of magnetic flux at high lati- 
tudes cause rotational histories which deviate up to 200% from 
those obtained following the WD approach. We find that this 
impact is large enough to mimic deviations of the dynamo ef- 
ficiency from linearity up to about 40%, and dynamo satura- 
tion limits at abou t 35 times the solar rotation rate (see also 
ISolankietal.il 9971) 

IWood & Linskvl Jl998l) and IWood et alJ d2002l) analysed 
the Lya absorption from the hot HI wall surrounding astro- 
spheres to determine mass loss rates of solar-like stars. The 
amount of absorption scales roughly with the square root of the 
ram pressure of the associated (magnetic) stellar wind, which 
is in turn a function of the mass loss and wind velocity. To con- 
vert the ram pressure into mass loss rates they assume a unique 
wind velocity (v w = 400 km/ s), taken to be isotropic and inde- 
pendent of the stellar rotation rate. Based on their results they 
determine an empirical relation between the mass lass rates of 
cool stars and their X-ray surface flux. In the framework of the 
present wind model the latitudinal profiles of the terminal wind 
velocity and ram pressure are however found to depend con- 



siderably on the stellar rotation rate. For slow rotators the ram 
pressure varies within a factor of ~ 3 with latitude. For rapidly 
rotating stars, in contrast, the latitudinal variation can rise to 
over two orders of magnitude. The latitudinal profile of the ram 
pressure in addition depends strongly on the non-uniformity of 
the surface magnetic field. If most of the open magnetic flux 
is located near the poles, then the strong magnetic fields there 
sustain the magneto-centrifugal driving at high latitudes despite 
the less efficient centrifugal forces close to the stellar axis of ro- 
tation; in this case the latitudinal variation of the ram pressure 
is about one order of magnitude. If most of the stellar mag- 
netic flux is concentrated at intermediate latitudes, then the po- 
lar magnetic field and wind acceleration will in contrast be very 
weak. Since the strong magneto-centrifugal driving at interme- 
diate latitudes entails high mass fluxes and wind velocities, the 
difference of the ram pressure with respect to the polar value is 
over two orders of magnitude. The wind ram pressure derived 
along a single direction is consequently not representative for 
the entire stellar surface. If rotational and magnetic particulari- 
ties of the star are not taken into account, the determined mass 
loss rates (and derived empirical relationships) can be subject 
to larger uncertainties. Since the difference of the latitudinal 
variation between a Polar Spot and a Latitudinal Belt distribu- 
tion is about one and half magnitudes, a considerable reduction 
of these uncertainties could be achieved by verifying which 
of the two cases describes the actual magnetic field distribu- 
tions of rapidly rotating stars more accurately. But for this, the 
yet small number of ZDI observations and field extrapolations 
has to be increased considerably. The X-ray luminosities of the 
stars considered by Wood et al. (2002) are below the X-ray sat- 
uration limit (Staufferet al. 1997) and the objects thus likely 
moderate rotators with rotation periods longer than about 2 d. 
The possible inaccuracies of the ram pressure should accord- 
ingly be less than an order of magnitude, and the uncertainty 
of the mass loss rates due to the latitudinal variation less then 
a factor of ~ 3. For stars of moderate rotation rate the linear 
empirical relationship between the stellar X-ray emission and 
mass loss rate should therefore hardly be affected. Due to the 
increasing latitudinal variation with rotation rate, the scatter of 
future observations of more rapidly rotating stars (with higher 
X-ray fluxes) around the predicted linear relationship is how- 
ever expected to increase. The conclusion of IWood et all that 
the mass loss rates of younger and more rapidly rotating stars is 
up to several orders of magnitude larger than in the case of the 
Sun, confirms the mass loss rates of rapid rotators determined 
in the framework of the WD approach. 

In their investigatio n of properties of stella r magnetic fluxes 
and magnetised winds, Schriiver et al 1 J2003I) combine obser- 
vational stellar data in a chain of arguments, which involves 
the functional dependence of the (effective) Alfvenic radius on 
the stellar rotation rate. Since in the framework of the present 
wind model, the difference between the effective Alfvenic ra- 
dius and its respective WD value hardly depends on the rotation 
rate, the impact of non-uniform field distributions on their work 
is expected to be marginal. In particular since their analyses is 
constrained to moderately active, slowly rotating stars, whose 
magnetic flux is typically located at low latitudes, in which case 
the effective Alfvenic radii are within ~ 10% comparable to the 
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WD value. However, as in more rapidly rotating stars magnetic 
flux is also found at intermediate and high latitudes, the tran- 
sition between different latitude-dependent field distributions 
implies an additional dependence of the Alfvenic radius on the 
rotation rate, which would have to be taken into account. The 
explicit functional dependence has to be determined observa- 
tionally (e.g. Strassmeie r 2002) or through numerical simula- 
tions (e.g. Sch ussler et alll996tlGranzer et all2000l) . 

In addition to the surface magnetic field distribution, the 
loss rates of magnetised winds depend significantly on coronal 
magnetic field topologies. Theoretical models indicate that the 
mass and AM loss rates of a multi-component corona consist- 
ing of both wind and dead zones (i.e., open and closed field 
regions, respectively) are considerably smaller than according 
to the WD model, since in addition to the confinement of coro- 
nal plasma in magnetic loops dead zones also alter the flow 
structure in adjacent w ind regions JPriest & Pneumar]ll974t 
iMestel & Spruidll987l) . Observations in the case of the Sun 
however indicate that the mass loss rate arising from the slow 
solar wind, traced back to closed equatorial field structures, is 
within a factor of two comparable with the mass loss rate of 
the fast solar wind, which originates from open coronal holes 
llWangll998l) . 

Simplified approaches to the wind model like the one 
used here typically ignore the non-linear influence of the 
trans-field component so that the resulting loss rates are sub- 
ject to inaccuracies. A comparison between in situ measure- 
ments of the solar wind by the Helios spacecrafts and re- 
spective WD-based estimates however showed that this ba- 
sic model is an adequate description of the physical interac- 
tion coupling the rotation of a star with its associated wind, 
insofar as simple magne tic stresses are ta ken as the princi- 
pal physical mechanism JPizzo et al.lll983l) : for a comparison 
between out-of-eclipse measurements by the Ulysses space- 
craft with a se lf-consistent model approach see, for exam- 
pic. ISautv el a l. (2005). Furthermore, investigations of the ro- 
tational evolution of stars, which make use of magnetic brak- 
ing timescales estimated according to WD-like wind mod- 
els, yield rotational histories which are consistent with ob- 
served distributions of stel lar rotation rates in young open clus- 
ters of different age ("e.g. lSolanki et aljll997l) . Whereas multi- 
dimensional MHD-simulations solve the trans-field component 
of the equation of motion self-con sistently (e.g., Sakurai 1985, 
Il990t iKennens & Goedbloedll2000l) . in the present work the 
poloidal component of the magnetic and velocity field are taken 
to be radial. The numerical simulations show that at large dis- 
tances from the star magnetic field lines do not remain radial 
but tend to be deflected away from the equatorial plane to 
higher latitudes, leadin g to a collimation of the wind a round 
the pole dTsineanos & Bogovalovll2000t IOkamotoll2000l) . The 
influence of this phenomenon on the mass and AM loss rates 
is not clear, since respective investigations are focused on the 
detailed wind structure of individual illustrative cases like the 
Sun. The combination of multi-dimensional MHD-simulations 
with more appropriate field distributions is desirable, although 
the large numerical demand makes such modelling more chal- 
lenging for stellar rotational evolution studies, where a large 
number of simulations is required to determine the relative im- 



pact of the evolving stellar structure, internal AM transport, dy- 
namo efficiency, dynamo saturation, and surface magnetic field 
distrib utions on the stellar rotation (e.g. iHolzwarth & Jardinel 
12004 . 

6. Conclusion 

The frequently observed non-uniform surface magnetic fields 
of active stars have a remarkable influence on stellar mass and 
AM loss rates. Concentrations of magnetic flux at high lati- 
tudes in the form of polar spots reduce the AM loss rate and ef- 
fective Alfvenic radius up to 60 % and 25% , respec tively, with 
respect to the approach of Web er & Davisl Jl967l) . The large 
deviations make them important ingredients for studies consid- 
ering the rotational evolution of stars, implying considerable 
consequences for the distribution of stellar rotation rates. The 
gathering of open magnetic flux in the form of a Latitudinal 
Belt has in contrast a strong impact on the error margins of 
empirical stellar mass loss rates derived using the ram pressure 
of stellar winds. In particular in the regime of rapidly rotating 
stars, the strong latitude-dependent magneto-centrifugal driv- 
ing causes latitudinal variations of the terminal wind velocity 
and ram pressure which span more than two orders of magni- 
tude. These examples show that the non-uniformity of surface 
magnetic fields represents an essential stellar property, whose 
neglect entails large uncertainties of observational and theoret- 
ical results. 

The loss rates and effective Alfvenic radii resulting from 
different non-uniform surface magnetic fields are found to 
cover a large range of values, although the total magnetic flux 
of each field distribution is the same. This shows that the clas- 
sification of stellar magnetic fields through a single, allegedly 
characteristic, field strength (e.g. the peak- or surface-averaged 
field strength) is considered to be inapplicable. Instead, possi- 
ble non-uniformities of the surface magnetic fields have to be 
taken into account, based either on empirical relationships or in 
the framework of theoretical constraints. In view of the impor- 
tance of the actual location of open magnetic flux, this makes 
an increase of the yet small number of combined Zeeman- 
DI observations and reconstructions of global stellar magnetic 
field topologies highly eligible. 
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Appendix A: Polytropic wind model 

The stellar wind is taken to be an ideal plasma with vanishing 
viscosity and infinite conductivity. Its structure is determined 
by the stationary equation of motion, and subject to a poly- 
tropic relation between the thermal pressure, p, and density, p. 
The polytropic index, F, quantifies the heating of the wind; a 
value close to one implies a nearly isothermal wind flow and 
a value close to the ratio of specific heats an almost adiabatic 
expansion. 

A.1. Basic relations and invariants 

In a reference system co-rotating with the stellar surface a sta- 
tionary plasma flow obeys the relation v co = KB/p, where v co is 
the flow velocity, B the magnetic field, and k the ratio between 
the mass flux and the magnetic flux; along a magnetic field line 
the latter is constant. The flow velocity in the rest frame of ref- 
erence is given by 



B 

v = v co + £2 x r = k h Q x r . 

P 



(A.l) 



where r is the vector pointing from the centre of the star to 
the location of a gas element, and Q a vector parallel to the 
stellar rotation axis, whose modulus, Q., is the stellar rotation 
rate. The wind is assumed to be axi-symmetric with respect to 
the rotation axis, en = e z , that is all derivatives in azimuthal 
direction, e$, vanish. 

Owing to the rotational symmetry, the azimuthal compo- 
nent of the equation of motion can be written in the form 



(A.2) 



with m — r sin 0. Substituting Eq. (IA.1> into Eq. (IA.2> yields 
that the specific angular momentum (per unit mass), 



1 

L = vatjj - - — Bj,m , 

4-7TK 



(A3) 



of plasma escaping along an open magnetic field line is con- 
stant. In addition to the common specific angular momentum 
of a moving gas element, oc r x v, this quantity also com- 
prises the tension of the magnetic field bent by the gas motion 1 . 
Following Eq. JA.3> . the azimuthal flow velocity is written in 
the form 



v 4> = 



(A.4) 



1 Considering the phy sical connection of Eq. <A.3> with the vortic- 
ity and current density, Goed bloed et"ail <2004l) use for the quantity 
L the nomenclature poloidal vorticity '-current density stream function 
and the symbol K to distinguish it from the specific angular momen- 
tum of gas flows in hydrodynamics. 
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For this expression to be finite at the Alfvenic point, (m A ,p A 
A-jik 2 ), the specific angular momentum has to be, 



L - Q.Tu 2 k = const. 



(A.5) 



The Alfvenic point marks the location along a magnetic field 
line where the Alfvenic Mach number is one (i.e. va = 
Z?a/ V^PaX an d the flow velocity changes from sub-Alfvenic 
to super- Alfvenic values. 

Since (in the rest frame of reference) the flow velocity of 
the wind is oblique to the magnetic field, the resulting Poynting 
flux leads to an energy transfer, 



— V [(VxB)xB] =£2v V(v cr) , 



(A.6) 



from the stellar rotation into the gas motion via the magnetic 
field, which contributes to the wind acceleration in the form of 
a magneto-centrifugal driving. Taking Eq. (IA.6> and the poly- 
tropic gas relation into account, the energy function, 



H = — + - 

2 r 



- 1 P 



(A.7) 



with m being the gravitational potential, has to be constant 
along a magnetic field line, 



H — E — const. 



(A.8) 



The four terms of Eq. iA.7\ are (in given order) the kinetic, the 
thermal, the gravitational, and the magneto-centrifugal contri- 
bution to the specific energy of the outflow. The Eqs. iA.H and 
JA.8> follow from the component of the equation of motion 
tangential to the wind flow and correspond to the Bernoulli in- 
tegral. The two constants, Eq. JA.5t and JA.8L are a direct con- 
sequence of the invariance of the problem against azimuthal 
and temporal translations (i.e. axial symmetric and stationary 
flow). For a m ore d etailed des cription of the theory of ma gne- 
tised winds see Mestel ( 1999): lLamers & Cassinellil (11999). 

A full treatment of the stellar wind problem requires the 
self-consistent solution of its trans-field component, which de- 
termines the poloidal field structure of the stellar wind. This 
aspect involves the treatment of a non-linear partial differential 
equation, which requires more sophistic ated theoretical and/or 
nume rical methods (e.g., Sakurai 1985; Vlahakis & TsinganosJ 
Il998l iKeppens & Goedbloedl l2000t) . More detailed investi- 
gations show that at large distances from the star, in the 
asymptoti c regime, magnetic field lines are deflected tow ard 
the pole JHevvaerts & Normanlll989HTsinganos & Bogovalovl 



2000; Okamoto 2000). In the present work the trans-field com- 
ponent is neglected and the poloidal field component instead a 
priori taken to be radial ("cf. lYehl 19761) . Consequently, the field 
lines form spirals around the rotation axis of the star, which 
are located on conical surfaces with constant openi ng angle. 
This a pproach is a simple generalisation of the lWeber & Davisl 
ill 9671) method in the sense that the individual boundary con- 
ditions and properties of open field lines are explicitly taken 
into account at all latitudes, although their mutual interaction is 
suppressed. Due to the latter simplification the absolute values 
should however only be considered as lowest-order estimates. 



A.2. Non-dimensional wind solutions 

The determina tion of wind solutions follows the method of 
ISakurail fl990); see also lSakurall l l 19851 Sect. 2). Using the non- 
dimensional variables 2 



y=£- , z=^-, (A.9) 

Pa vV,a 



Eq. iA.li is written in dimensionless form, 



r A Kb , K n ( xzU-X 1 ) 

H = H—Z- = —z 2 + — — 

GM, 2 2 1 - x 2 z 
v 

1 K T , 2 u-r 
with the non-dimensional parameters 



r l l xU 

Kt = T. f~T ' Ka = — ' and K P = 



PaJo 



A)*Pa 



(A. 10) 



(A.ll) 



Quantities with index '0' refer to values at the reference level, 
ro- 

The parameter Kq is the only quantity that mediates a de- 
pendence of the wind on both the stellar rotation rate and lat- 
itude. In the equatorial plane the magneto-centrifugal driving 
of the stellar wind is quantified through the co-rotation radius, 
r co = -\/GM t /Q. 2 , which specifies the distance from the star 
where the inward directed gravitation is balanced by the out- 
ward directed centrifugal force. Here, the magnetic field lines 
outside the equatorial plane are situated on coni with constant 
opening angles. Hence the latitude- and rotation-dependent ra- 
dius 



CM, 

a 2 



sin 



-2/3, 



(A. 12) 



denotes the point beyond which in radial direction, e r , the cen- 
trifugal acceleration of the gas, a c , outbalances the gravita- 
tional acceleration, g, that is (g + a c ) • e r = 0. The param- 
eters Kj and depend on the boundary conditions xo and 
yo = 1/(XqZq), determined at the reference radius, ro. There, 
the base temperature, Tq, of the wind defines the Parker radius, 



'"Pa 



GM t GM*f£ 1 



2c 2 



221 T 



(A.13) 



where c^o is the (isothermal) sound speed, and ji the mean 
atomic weight of the plasma. The dependence of the wind on 
the magnetic field strength enters through the parameter 



forgo 
Bl 



(A. 14) 



which denotes the ratio between the gas pressure and the mag- 
netic pressure at the base of the wind. 



2 Since along a magnetic field line the mass flux per solid angle, 
- r 2 pv r , is constant, the relation ^yz = 1 holds. In contrast to 



Sakurai Il985h . here the flow velocity is used as an independent vari 
able instead of the density. 
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Since 'H = E — const., a non-dimensional wind solution, 
z(x), is represented in the (x, z)-plane by an isoline. Along this 
trajectory the total differential of the energy function, 



"-(£).*♦(£ 



dz ■ 



(A. 15) 



vanishes, so that the variation of the flow velocity in radial di- 
rection is determined by the differential equation 



dz 
dx 



m_ 

'(f)/ 



(A. 16) 



To prevent this derivative from becoming infinite at points 
where the numerator vanishes, the denominator simultaneously 
has to be zero. The condition 



(A.17) 



defines two locations along the trajectory, which are usually re- 
ferred to as the slow and fast critical points, (x s , z s ) and (xf, zi), 
respectively. Since both points are situated on the same energy 
level, the condition 



(A.18) 



holds. The five non-linear algebraic Eqs. (IA.17I for the slow 
and fast critical points) and JA.18> are used to determine sets 
of non-dimensional wind solutions (x s ,z s , Xf,Zi, Kp) as a func- 
tion of the two remaining parameters Kj and Kn. These solu- 
tions are as yet independent of specific boundary conditions; 
the only physical quantity explicitly entering Eq. 



JA.10> is the 

polytropic index. For a given value of T, a look-up table is pre- 
calculated to speed up the following determination of dimen- 
sional wind solutions, which depend on the boundary condi- 
tions at the reference level. 



choice of Kj and K&\ note that Kp and £ are both implicit func- 
tions of these independent variables. Once a non-dimensional 
wind solution has been found, the Alfvenic radius and density 
are determined by back-substitution of the parameters Kj, 
and Kp. 

A.4. Terminal wind velocity 



The terminal wind velocity, Zoo, is determined using Eq. JA.10> 
in the limit of large distances from the star (x — > oo), which 
yields the cubic equation 



z 3 -z2 



K 8 



1 + A +2^ = 0. 
Kn) Kp 



(A.21) 



The correct solution is subject to the constraint that even far 
away from the star the wind velocity has still to increase, that 
is 

2z 2 K,y Y - 1 



dz 
dx 



x {z'Kp-Kn) 



>0. 



(A.22) 



For a stationary outflow to exist, the terminal wind velocity has 
to exceed a latitude-dependent minimal flow velocity, 



> z M 



r^Q VM . 2/3 n 
I = sin /J 6 . 

Vr,A 



(A.23) 



with vm = -^ Q 2 r^ v n A being the characteristic Michel veloc- 
ity of the wind (Michel 1969). With a proper set of parameters 
Kn,Kp, and E obtained from Eqs. dA. 19i and JA.20> . the ter- 
minal velocity is determined from Eqs. (IA.2H and (IA.23> . 



A. 3. Dimensional boundary conditions 

In addition to the stellar rotation rate, Q., a dimensional wind 
solution requires the definition of the temperature, 7q, the den- 
sity, po, and the radial magnetic field strength, B n o at the ref- 
erence level, ro. The link between these dimensional boundary 
conditions and the dimensionless wind solutions above is es- 
tablished through a combination of Eqs. JA.11> . 

*4-«^.!(£P(a)V. <A.i9> 

to eliminate the yet unknown Alfvenic radius, r\, and density, 
Pa, as well as the energy relation 

The latter is determined at the reference level, where the spe- 
cific kinetic energy of the wind flow is assumed to be negligible 
compared with the thermal and centrifugal energy. For a given 
set of boundary conditions, the right sides of Eqs. dA. 191 and 
JA.20> are constant and the solution found by an appropriate 



